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Abstract 

In our paper Phys. Rev. Lett. 92, 151801 (2004), the oscillations of 
strongly deformed Bogomolny-Prasad-Sommerfield magnetic monopoles 
have been studied on a fixed Minkowski background. The purpose of 
the present article is to provide a more detailed account on the results 
yielded by our numerical simulations. In particular, an analysis on the 
dependence on the strength of the initial excitation is carried out in order 
to distinguish features which are already present for small perturbations 
from those which are consequences of the nonlinearity of the system. 



1 Introduction 

The investigation of solitons in particle physics is of fundamental interest 
(see e.g. for a recent review). In particular, considerable attention has 
been paid to the study of 't Hooft-Polyakov magnetic monopole solutions 
of coupled Yang-Mills-Higgs (YMH) systems |5J|3|. Magnetic monopoles 
are important because they are present in a wide range of gauge theory 
models. Moreover, their existence can be used to provide an account on 
various physical phenomena, including e.g. the quantized nature of electric 
charge. In a recent paper by the present authors 0] a report is provided 
on the numerical study on the behavior of strongly excited Bogomolny- 
Prasad-Sommerfield (BPS) magnetic monopoles 5 6 . In a related article 
by Forgacs and Volkov |7j the same system is investigated by perturbative 
methods. There is a striking consistency between the key results of these 
two independent studies. The aim of the present paper is to provide more 
detailed description of the system based on the results obtained by the 
numerical investigations of 0]. 

The investigated dynamical magnetic monopole is a solution of a cou- 
pled SU(2) YMH system. The Yang-Mills field is represented by an su(2)- 
valued vector potential A a and the associated 2-form field F a b reads as 
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Fab = VaA b - VbAa + ig [Aa, Ab] (1) 

where [ , ] denotes the product in su(2) and g stands for the gauge cou- 
pling constant. The Higgs field (in the adjoint representation) is given by 
an su(2)-valued function ip while its gauge covariant derivative reads as 
V a ip = V a tp + ig[A a ,ip]. The dynamics of the investigated YMH system 
is determined by the Lagrangian 

£ = Tr(F ef F ef ) + 2Tr(V e ipV e iP) + A [TrftWO - v 2 } 2 , (2) 

where A is the self interaction constant of the Higgs field. In order to have 
finite energy solutions the value of Tr(ipip) must tend to v 2 as the distance 
from the monopole goes to infinity. The energy-momentum tensor of the 
YMH field takes the form 



T - 1 
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Tr(FaeF b e ) - Tr{Vai>V b ^) + jg ab £ 



(3) 



Our considerations are restricted to spherically symmetric configu- 
rations yielded by the 'minimal' dynamical generalization of the static 
't Hooft-Polyakov magnetic monopole configurations [2111] (see also [8]). 
Accordingly, the evolution takes place on Minkowski spacetime the line 
element of which, in spherical coordinates (t, r, 6, </)), is 

r 2 - r 2 (d6 2 + sin 2 9 dcj> 2 ) , (4) 

while the Yang-Mills and Higgs fields, in the so called abelian gauge, are 
assumed to posses the form 

Aa = -- [w{r 2 (d6) a - t x sin 0(dct>) a } + T 3 cos 9{d(j))a] (5) 

i> = Hr 3 , (6) 

where the generators {tj} (1=1,2,3) of su(2) are related to the Pauli ma- 
trices a 1 as Tj = hffj, moreover, w and H are assumed to be smooth 
functions of t and r. 

In the Bogomolny-Prasad-Sommerfield (BPS) limit, i. e. when A — > 0, 
the mass of the Higgs boson becomes much less than the mass of the 
vector boson. The field equations in the A = case become 

r 2 d 2 w — r 2 d 2 w = w [(w 2 — l) + g 2 r 2 H 2 ] (7) 
r 2 d' 2 H + 2rd r H - r 2 d 2 H = 2w 2 H , (8) 

with the restriction that 

lim H = v , (9) 

r — >oo 

in order to represent the BPS limit of finite energy solutions. We note 
that the time independence of the limit value of H at both spacelike and 
null infinity follows from the field equations |7) and © (see |0])- 

The system of equations and © with condition JUJ has an analytic 
solution, the static BPS monopole |5J|Sj 



w o = -^—n r > Ho = v 

smh(gvr) 
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In the rest of the paper considerations will be restricted to the A = 
case and strong impulse type excitations of the BPS monopole will be 
investigated. 

In the BPS limit the Higgs field becomes massless and the only scale 
parameter of the system is the vector boson mass m w = gv. Since in the 
case considered here v 7^ 0, the rescalings i — > t = tm w , r — > f = rm m and 
H — > H — H/v transform the parameters to the value g = v — 1. This 
implies that it suffices to consider the evolution of the g = v = 1 system 
numerically and to study merely the dependence of the evolution on the 
various initial conditions for this system. 

To have a computational grid covering the full physical spacetime - 
ensuring thereby that the outer grid boundary will not have an effect on 
the time evolution - the technique of conformal compactification, along 
with the hyperboloidal initial value problem, is used. This way it is pos- 
sible to study the asymptotic behavior of the fields close to future null 
infinity, as well as, the inner region for considerably long physical time 
intervals. 

The conformal transformation we use is a slight modification of the 
static hyperboloidal conformal transformation applied by Moncrief |10|. 
It is defined by introducing first the new coordinates T and R instead of 
t and r as 

T = ujt+l~ ^/iu 2 r 2 + 1 and R= ^ 2r2 + J__1 (n) 

ujr 

where iv is an arbitrary positive constant. The Minkowski spacetime is 
covered by the coordinate domain given by the inequalities — 00 < T < 
+00 and < R < 1. The R = const, lines (same as the r = const. 
lines) represent world-lines of 'static observers', while the T = const. 
hypersurfaces are hyperboloids of the Minkowski spacetime satisfying the 
relation (cut + 1 — T) 2 — uj 2 r 2 — 1. The constant u is introduced here in 
order to have control on the size of the monopole region in terms of the 
R coordinate, which is useful for the optimization of the numerical code. 
In the rest of the paper we present results obtained by setting ui = 0.05. 

The line element of the conformally rescaled metric g a t = Q 2 g a b in the 
coordinates (T, R, 8, (f>) takes the form 

Q 2 

dt = —dT 2 + 2RdTdR - dR 2 - R 2 {d6 2 + sin 2 6» d(f> 2 ) (12) 

where the conformal factor is 

Q, = u(l - R 2 )/2 . (13) 

The R = 1 coordinate line represents through which the conformally 
rescaled metric g a b smoothly extends to the unphysical coordinate domain 
with R > 1. The simple relation rQ. — R holds between the old and new 
radial coordinates. 

Using the substitution H(t,r) = h(t,r)/r + v the field equations J7J 
and © in the new coordinates read as 

q3w = ™[(w 2 -l) + g 2 (h + vRQ- 1 ) 2 ^ (14) 

«p/i = 2 (h + vRfl^ 1 ) w 2 , (15) 
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where the differential operator <p is defined as 



4R 

V = 



(R 2 + l) 2 



— ^8%, - &t - 2Rd R d T 

LU 2 



2tt VtR{R 2 + 3) 



uj{R 2 + \) uj(R 2 + 1) 



(16) 



These equations can be put into the form of a first order strongly hyper- 
bolic system 0. The initial value problem for such a system is known to 
be well-posed |11|. In particular, we solved this first order system numeri- 
cally by making use of the 'method of line' in a fourth order Runge-Kutta 
scheme following the recipes proposed by Gustafsson et al All the 

details related to the numerical approach, including representations of 
derivatives, treatment of the grid boundaries are to be published in 
The convergence tests justified that our code provides a fourth order rep- 
resentation of the selected evolution equations. Moreover, the monitoring 
of the energy conservation and the preservation of the constraint equa- 
tions, along with the coincidence between the field values which can be 
deduced by making use of the Green's function and by the adaptation of 
our numerical code to the case of massive Klein-Gordon fields, made it 
apparent that the phenomena described below have to be, in fact, physical 
properties of the magnetic monopoles. 

In each of the numerical simulations initial data on the T = hypersur- 
face was specified for the system of our first order evolution equations. In 
particular, a superposition of the data associated with the BPS monopole, 
1101 . and of an additional pulse of the form 

-exp [t— fz^p-l : i£ r S[a-b,a + b]; 
Uhtr) - { (17) 

, otherwise, 

with a > b > 0, which is a smooth function of compact support, was 
used. This choice, providing non-zero time derivative for w, corresponds 
to "hitting" the static monopole configuration between two concentric 
shells at r = a — b and r — a + b with a bell shape distribution. Basically 
the same type of evolution occurs when instead of (dTw)o we prescribe 
(9rft)o in a similar fashion. The only distinction is that in this later case 
the energy of the pulse is given directly to the massless Higgs field which 
carries immediately a considerably large portion of it to ,J? + before the 
Yang-Mills field could take it over to feed it into the breathing state of the 
monopole and into the expanding oscillations of the massive Yang-Mills 
field. 

All the simulations shown below refer to the same pulse shape 1171 
corresponding to the choice of the parameters a — 2, b — 1.5 and d = 10, 
only the amplitude c is varied. It is important to keep in mind that the 
energy of the pulse is proportional to c 2 . For example, for the case c = 70 
presented in |lj it is 55.4% of the energy of the static monopole, which is 
clearly not just a simple perturbation. We would also like to emphasize 
that the figures shown below are typical in the sense that for a wide range 
of the parameters characterizing the exciting pulse, qualitatively, and in 
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certain cases even quantitatively, the same type of responses are produced 
by the monopole 

In order to provide a better picture about the size and distribution of 
the monopole and its relation to the initial data it is instructive to plot 
the energy density distribution on the T — initial hypersurface. The 
conserved energy current vector j a = T a b u a is formed from the energy 
momentum tensor T a b , where u = is the velocity vector of the the static 
observers characterized by the constant value of r along their worldline. 
The projection s = j a n a = luT° (R 2 + l)fi _1 /2, where n a is the future 
pointed unit normal to the constant T hypersurfaces, gives a gauge in- 
variant quantity characterizing the energy density on these surfaces. The 
energy "density" associated to shells of radius R can be calculated as 
£ = J E^/\h\d6d<j) = 2ttujT° q (R 2 + l)R 2 n~ 4 , where h is the determinant 
of the induced metric on the constant T hypersurfaces. On Fig.Qwe plot 
the value of e and £ for the static monopole both with and without an 
initial deformation of the form of 1171 . The plot of £ is more instructive 
because its integral from the center R — to null infinity R = 1 (i. e. 
the area below the curve) gives the total energy of the system along the 
chosen constant T hypersurface. As the time passes, the value of this 
integral decreases exactly by the energy radiated to null infinity. 




Figure 1 : Energy density distributions e and £ on the T = initial hypersurface 
for the static BPS monopole (c = 0) and for two different choices of the initial 
amplitude c. 



Choosing temporarily c = 70, we consider next the quantity w, corre- 
sponding to the massive Yang-Mills field, plotted on succeeding T = const. 
hypersurfaces, providing thereby a spacetime picture of its time evolution 
(see Fig. |2l . One can see the formation of expanding shells of high fre- 
quency oscillations in the asymptotic region |12) . These shells take out 
energy from the central monopole region, although they never reach null 
infinity because of their massive character. The time evolution of the 
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Figure 2: Spacetime diagram showing the time evolution of the variable w 
corresponding to the massive Yang-Mills component of the gauge field for the 
time interval < T < 3.71. The amplitude of the oscillations can be evaluated 
by observing that at the center w = 1, while at infinity w = 0, always holds. 



other gauge variable, the massless Higgs field h is strikingly different (see 
Fig-BJ- There is no sign of the high frequency oscillations in h, justifying 
that the two fields effectively decouple in the asymptotic region. Since the 
oscillations in h are nonzero at null infinity, the direct energy transport 
to J? + by the Higgs field, with the velocity of light, is apparent. 

From Fig.|Uand Fig-Hit appears that the lower frequency oscillations 
are not concentrated only at the monopole but they are present basically 
everywhere. However, if one considers the clearly physically meaningful 
quantity, the energy density £ corresponding to spherical shells of radius 
R, it is apparent that the oscillations are concentrated at the central 
region, forming a long lasting 'breathing state' of the monopole. In Fig. 
2] we show the energy density difference £ — £ o in a large central region, 
where £o is the energy density of the static monopole. (See Fig. 1 of 
for the plot of £ up to null infinity.) Although Fig. 2] covers a relatively 
early time period, the picture remains essentially the same for the rest of 
the time evolution when taking time intervals of the same length, with 
the only difference of an overall decrease of its amplitude proportional to 
T~ 5/6 (see Fig. 5 of |U). 

Next we consider the dependence of the energy density oscillations on 
the amplitude of the initial deformation. The energy density p — T a bU a u 
observed by static Minkowski observers cannot be used to form a con- 
served quantity when integrating on the hyperbolic T — const, hypersur- 
faces. However, for our specific choice u = 0.05 the T = const, hyper- 
boloids are so close to the flat t = const, surfaces in the inner monopole 
region that plots of the quantities e and p are almost indistinguishable. 
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Figure 3: Spacetime diagram showing the time evolution of the variable h cor- 
responding to the massless Higgs component of the gauge field for the initial 
amplitude choice c = 70. At the center h = for all time, while at infinity h 
remains —1 until the radiation arrives there. 
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Figure 4: Spacetime diagram showing the time evolution of the energy density 
difference £ - S horn T = 4.97 to T = 5.58. 



Furthermore, for the static BPS solution e = p = po. On Fig. Q>]we plot 
the energy density fluctuation p — po around the static monopole value po 
for four different choices of the amplitude c. Since for small initial per- 
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turbations the deviations of the density p from the static value po appear 
to be proportional to the amplitude of the perturbation, we rescale the 
amplitude of the four graphs accordingly, by plotting the value (p — po) /c. 
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Figure 5: Evolution of the rescaled energy density difference (p — po)/c as seen 
by a constant radius observer at R — 0.0254 (i. e. r — 1.018) for four different 
values of c and two different time intervals. 

For c = 0.7 and c — 7 the corresponding curves almost entirely overlap, 
showing that these excitations are still in the linear domain. Indeed, the 
energy provided by the initial perturbation in those two cases are 0.00554 
and 0.554 percent, respectively, of the energy 4-7T of the static monopole. 
For the c — 70 case the energy given to the system is 6.966, which is 
55.4 percent of the energy of the initially static monopole. In this case 
the behavior is clearly nonlinear initially, but after the direct pulses have 
left the system the quasinormal oscillation of the monopole is still not 
too far from the linear regime. For c = 280 the additional energy is 
111.46, which is more than 8 times the energy of the static monopole. In 
this case even the monopole oscillations show a different character at the 
beginning, having a significantly lower frequency. However, in all cases, 
the frequency of the monopole oscillations tend as t~ 2 / 3 to the asymptotic 
value 1 corresponding to the vector boson mass in our rescaled coordinates 
(see Fig. 5 of g|). 

After the initial period, i.e. for t > 5, the monopole starts a quasi- 
periodic oscillation with amplitude decreasing as t~ 5 ^ 6 . It is visible that 
the mean value of these oscillations, which represents the average energy 
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density at the given radius is negative for large initial deformations. This 
means that the energy contained in the oscillating monopole region is ac- 
tually smaller than the energy of the static monopole which is a clear man- 
ifestation of the nonlinear character of the evolution. As it was pointed 
out to us by Michael Volkov, this situation is similar to that where a large 
stone is dropped into a lake, and the water level decreases because of the 
submerging of the stone. This analogy, however, does not provide a com- 
pletely satisfactory explanation in the monopole case, because the average 
energy remains negative during the entire evolution, approaching merely 
in the asymptotic limit the static value from below. This might lead to 
the interpretation that the monopole behaves like an extremely viscosal 
fluid. 

In the last part of the paper we consider the energy radiated to null 
infinity up to the time T. This energy can be calculated as E r = f SdT 
where S = lim, — >00 j a k a with k a being the unit normal to the constant 
r hypersurfaces. The possibility that we could study this radiation, in 
particular the evolution near J? + , is due to the fact that the hyperboloidal 
initial value problem was used in our numerical simulations. The radiation 
arrives to when the light cone emanating from the outer edge of the 
initial perturbation reaches null infinity, i.e. approximately at T = 0.86. 
For small perturbations the radiated energy is proportional to the square 
of the initial amplitude. Hence we plot the rescaled energy E r /c 2 on Fig|S| 
for various initial data specifications. 

The energy given to the system by the initial deformation is approx- 
imately 0.0014216 • c 2 . For small initial perturbations the total radiated 
energy is 0.00068273 • c 2 , which is 48.02 percent of the energy given to 
the system. The rest of the energy is in the self-similarly expanding shells 
formed by the massive Yang-Mills field w. These shells get arbitrarily far 
from the monopole region, but because of their massive character their 
velocity is smaller than the speed of light, whence they never reach null 
infinity. Actually, the energy in the shells is initially a little bit larger 
than the remaining 51.98 percent, since as we have seen, the average en- 
ergy in the central monopole region is smaller than the energy of the static 
monopole solution in the same region. As the time passes the amplitude 
of the monopole oscillations decreases and the mean value of the energy 
density approaches the static value from below. Accordingly, there is, al- 
though a low scale, energy transfer back to the monopole from the outer 
regions ensuring, in turn, that in the inner region the system will settle 
down to the static monopole solution. 

We would also like to emphasize that - as it is intuitively expected - 
the ratio between the energy carried by the shells and the energy radiated 
to J? + depends on the profile of the initial excitation. It is found that in 
case of the specific choice 1171 this ratio depends on a, b and d, but for 
small amplitudes in the linear regime it is apparently independent of the 
parameter c. 

Various new features of the evolution of excited 517 (2) BPS monopoles 
have been presented here. The results contained in the last part of this 
paper indicate that there remained a lot of interesting issues to be stud- 
ied both numerically and analytically in the dynamical behavior of the 
t'Hooft-Polyakov monopoles. 
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Figure 6: The rescaled energy E r /c 2 radiated by the system until time T for 
various initial amplitudes c. On the lower figure the logarithmic plot of the 
difference from the asymptotic value Ef is displayed for a longer time period, 
indicating that the energy tends to Ef as T~ 2 / 3 . 



Acknowledgments The authors wish to thank Peter Forgacs who 
suggested the numerical investigation of magnetic monopoles. We also 
would like to thank him, Michael Volkov and J6rg Frauendiener for dis- 
cussions. This research was supported in parts by OTKA grant T034337, 
TS044665 and NATO grant PST.CLG. 978726. 



References 

[1] V. Rubakov: Classical Theory of Gauge Fields, Princeton Univ. Press 
(2002) 

[2] G. 't Hooft: Magnetic monopoles in unified gauge theories, Nucl. Phys. 
B 79, 267-284 (1974) 

[3] A.M. Polyakov: Spectrum of particles in quantum field theory, JEPT 
Lett. 20, 194-195 (1974) 



Fodor and Racz 



[4] G. Fodor and I. Racz: What Does a Strongly Excited 't Hooft-Polyakov 
Magnetic Mono-pole Do? Phys. Rev. Letters 92, 151801 (2004) 

[5] E.B. Bogomolny: Stability of classical Solutions Sov. J. Nucl. Phys. 
24, 449 (1976); Yad. Fiz. 24, 861-870 (1976) 

[6] M.K. Prasad and CM. Sommerfield: Exact classical solution for the 't 
Hooft Monopole and the Julia-Zee dyon, Phys. Rev. Lett. 35, 760-762 
(1975) 

[7] P. Forgacs and M.S. Volkov: Resonant Excitations of the 't Hooft- 
Polyakov Monopole Phys. Rev. Letters 92, 151802 (2004) 

[8] K. Huang: Quarks, leptons & gauge Gelds, World Scientific Pub., 
Singapore (1992) 

[9] G. Fodor and I. Racz: Numerical investigation of the time evolution 
of Yang-Mills-Higgs systems, in preparation 

[10] V. Moncrief: Conformally regular ADM evolution equations, talk 
given at the workshop on Mathematical Issues in Numerical Relativity 
held at ITP Santa Barbara from January 10-28, 2000, for an online 
version see http : / /online .ftp .ucsb .edu /online /numrelOO /moncrief/ 

[11] B. Gustafsson, H.-O. Kreiss and J. Oliger: Time dependent problems 
and difference methods, (Pure and Applied Mathematics) New York: 
Wiley (1995) 

[12] G. Fodor and I. Racz: Massive fields tend to form highly oscillating 
self-similarly expanding shells , Phys. Rev. D. 68, 044022 (2003) 



